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ABSTRACT 

We employ a bias-corrected abundance matching technique to investigate the coevo- 
lution of the ACDM dark halo mass function (HMF), the observationally derived 
velocity dispersion and stellar mass functions (VDF, SMF) of galaxies between z — 1 
and 0. We use for the first time the evolution of the VDF constrained through strong 
lensing statistics by Chae (2010) for galaxy-halo abundance matching studies. As a 
local benchmark we use a couple of z ^ VDFs (a Monte-Carlo realised VDF based 
on SDSS DR5 and a directly measured VDF based on SDSS DR6). We then focus on 
connecting the VDF evolution to the HMF evolution predicted by A^-body simulations 
and the SMF evolution constrained by galaxy surveys. On the VDF-HMF connection, 
we find that the local dark halo virial mass-central stellar velocity dispersion (Mvir-c) 
relation is in good agreement with the individual properties of well-studied low-redshift 
dark haloes, and the VDF evolution closely parallels the HMF evolution meaning lit- 
tle evolution in the M^ir-cr relation. On the VDF-SMF connection, it is also likely 
that the stellar mass-stellar velocity dispersion (Af^-cr) relation evolves little taking 
the abundance matching results together with other independent observational results 
and hydrodynamic simulation results. Our results support the simple picture that as 
the halo grows hierarchically, the stellar mass and the central stellar velocity disper- 
sion grow in parallel. We discuss possible implications of this parallel coevolution for 
galaxy formation and evolution under the ACDM paradigm. 

Key words: galaxies: evolution - galaxies: formation - galaxies: haloes - galaxies: 
kinematics and dynamics - galaxies: statistics - galaxies: structure 



1 INTRODUCTION 



The current ACDM hierarchical structure formation the- 
ory predicts robustly the evolution of the dark halo mas^ 
funct i on (HMF) ov er cos mic time ( e .g., [S pringcl ct al. 

l2007l : iLukic et al.l 
2010l ). Because vis- 



2005; Warren et a l. 2006: Rccd ct al 
2007. : Tinker et ah ,200& : .Kl ypin et al"^ 



ible galaxies are believed to form and reside within the 
haloes under the ACDM paradigm, the statistical func- 
tions of galaxies, such as the luminosity function (LF), the 
stellar mass function (SMF), and the stellar velocity (dis- 
persion) function (VF, VDF), are also expected to evolve. 
Connection of these statistical functions of galaxies with 
the theoretical HMF is not straightforward due to, and 



^ Throughout this refers to the total mass within the virial radius 
of the halo. Accordingly, it includes the stellar mass once the 
galaxy is formed. 



mirrors, the complex processes of galaxy formation and 
evolution including star formations, supernovae explosions, 
AGN activities and galaxy merging. Some recent works in 
the literature are focused on the connection of the HMF 
with the broa dly measured stellar mass function (SMF) of 
galaxies (e.g.. IConroy fc Wechsleij|2"009l : iMoster et ai]|2009l : 
IGuo et al.ll2009l ). Notice that the SMF as well as the LF have 
mainly to do with the star formation history of galaxies. 

Galaxy formation in the halo has dynamical conse- 
quences as well. As stars are formed in the inner halo, the 
halo responds and the inner halo dark matter distribution is 



modified (e . g.. iBlumenthal et al.l 119861: iGnedin et al 112004 
iRudd et all l2008l : lAbadi et al.1 bOOalTissera et al.1 boiol ). 



Consequently, not only the total (i.e. dark plus stellar) mass 
distribution but also the dark matter distribution may be- 
come different from the pure dark matter distribution pre- 
dicted by the ACDM. This dynamical aspect of galaxy for- 
mation is a crucial part of cosmological studies. Ultimately, 
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the theory of galaxy formation should predict successfully 
the dynamical evolution as well as the star formation history 
of galaxies. The statistical property of the dynamics of galax- 
ies is encoded in the VDF of galax i es. Th e local total VDF 
is carefully reconstructed bv lChad (|2010l ') using SDSS DR5 
galaxy counts and intrinsic correlations bet ween luminosi- 
ties an d velocities of galaxies. More recently. iBernardi et al.l 
estimates directly the local total VDF based on SDSS 
DR6 measurements (of a DR4 sample). Chaa ((201Q ') then 
constrains the evolution of the VDF up to 2 ~ 1 through 
the statistical properties of strong lensing galaxies based on 
the empirical result that the average total (luminous plus 
dark) mass profile of galaxies is isothermal in the optical 
region. IChael(|2010l ) notices that the differential evolution of 
the derived VDF is qualitatively similar to the evolution of 
the theoretical HMF under the current ACDM paradigm. 

In this work we make a detailed quantitative compar- 
ison between the VDF evolution constr ained f rom s trong 
lensing statistics through the method of IChad l|201(t ) and 
the evolutions of mass functions (i.e. the HMF and the SMF) 
from the literature. In doing so, we investigate the local (sta- 
tistical) correlations of the velocity dispersion (a) of a galaxy 
with the virial mass (Mvii) of the surrounding halo and the 
stellar mass (M*) of the galaxy, i.e. (j(Mvir) and (t(M*), 
and their evolutions out to 2 ~ 1. These empirical correla- 
tions will provide independent statistical constraints on the 
structures of galaxies and haloes and their evolutions. We 
investigate the implications of the local correlations for the 
baryon-modified dark halo structures in a following work. 
In this work we focus on the evolutions of the correlations. 
We find that the halo mass, the stellar mass and the stel- 
lar velocity dispersion are coevolving in a parallel way for 
^ 2 < 1. We discuss the implications of this result for 
galaxy formation and evolution under the ACDM paradigm. 

This paper is organised as follows. In 2, we describe 
the method of analysis and the statistical functions to be 
used in this work; some details of the analysis are given in 
the Appendix A. In §3, we investigate the connection be- 
tween the (evolving) VDF of galaxies from the SDSS and 
strong lensing statistics and the HMF from N-body simu- 
lations. We obtain the relation cr(Mvir) and its evolution. 
We also examine the compatibility of the evolutions of the 
VDF and the HMF. In §4, we compare the VDF evolution 
with the SMF evolution from galaxy surveys. We investigate 
the evolution of a{A'h) and the compatibility of the cur- 
rent observationally constrained VDF and SMF. In §5, we 
discuss the implications of the results for galaxy formation 
and evolution in the context of the current ACDM struc- 
ture formation paradigm and cosmological observations. We 
conclude in §6. Unless specified otherwise, we assu me a 
WMAP 5 year ACDM cosmology (|Dunklev et al.ll2009l ) with 
(fimcf^Ao) = (0.25,0.75) and Ho = lOO/ikms"! Mpc"\ 
When parameter h does not appear explicitly, /i = 0.7 is as- 
sumed. In Appendix B we compare t he results of th i s wor k 
with the near l y con current results bv lDutton et al.l l|2010l ). 
iDutton et all i|2010l ) focus on the connection between the 
circular velocity in the optical region (at about the pro- 
jected half -light radius) Uppt a nd that at the virial radius 
Vvir- While [Putton et al.l (|2010l ) use various estimates of the 
halo mass including satellite kinematics, weak lensing and 
abundance matching, their results ar e confined to 2: ^ . 
Another important difference between iDutton et al.l (|2010l ) 



and this work is that button et all (|2010l ) use observed stel- 
lar mass- velocity relations while this work uses the observa- 
tionally derived velocity dispersion functions for abundance 
matching. 



2 METHOD: ABUNDANCE MATCHING OF 
STATISTICAL FUNCTIONS 

A statistical function of galaxies or haloes at a given epoch 
is defined by 



(x) = 



dn{> x) 



dx 



(1) 



where x is the variable under consideration (e.g. a, M*, 
Mvir) and n(> x) is the integrated comoving number density 
down to X. 

We use the a bundance matching i nethod (e.g. 
Kravtsov et al. I I2OO4I : IVale fc Ostrikeil |2004| : IConrov et al.l 



2006h to relate statistically one variable {x) to another {y). 



x{y) from n(> x) = n(> y). 



(2) 



Namely, we have 

y — y{x) or x - 

The key assumption for equation ([2} to be valid is that the 
two variables are monotonically increasing functions of each 
other. The accuracy of the median relation between x and y 
derived from the abundance matching method depends on 
the nature of the intrinsic scatter of the t rue relation (see 
iTasitsiomi et al.ll20o3 : lBehroozi et al.ll2010l '). In Appendix A, 
a simulation is carried out to investigate the possible effect of 
the intrinsic scatter. It turns out that based on an observa- 
tionally motivated intrinsic scatter the abundance matching 
method reproduces the overall behaviour of the intrinsic re- 
lation up to a maximum bias of ~ 0.08 dex in most cases. We 
estimate and correct the biases in our abundance matching 
analyses. 

This work is primarily concerned with connecting the 
stellar velocity dispersion of a galaxy (a) to the dark halo 
virial mass (Aivir) and the galaxy stellar ma ss (M <. ). For 
the VDF at 2 = we u se the resuhs from IChael (I2010l ) 
and IBernardi et all 1 2010 (Fig. [T]). Specifically, we use the 
'AO' VDF of IChael (120101 1 that is a resuh firom combining 
the early-t ype and the lat e -type VDFs based on SDSS DR5 
data. The IBernardi et all ([2010') VDF is a direct fit to all - 
type galaxies based on SDSS DR6. Bernardi et al.l (|2010| ) 
give various fit results depending on the range of velocity 
dispersion and fit method. For this work we use the fit re- 
sult for a > 125 km s~^ (taking into account measurement 
errors) because this work is primarily concerned with the 
evolution of massive galaxies. W e then use strong lensing 
statistical analysis of lChad (|2010l ) to constrain the evolution 
of the VDF u p to z = 1. T he results are shown in Fig. [1] The 
resuh f or thelchad (I2OI0I) VDF is a reproduction while that 
for the IBernardi et all (120101') VDF is a new result. Notice 



that for the Bernardi et al.l (l2010t) VDF a m odified Schechter 
fun ction i n trodu ced by Sheth et al.l (|2003h is used while for 
the IChad l|2010t ) VDF a correction term is included. The 
simplicity of the IBernardi et al] l|2010l ) VDF allows all four 
parameters of the function to be varied and constrai ned by 
strong lensing data in contrast to the IChad (|2010D VDF 
for w hich some parameters must be fixed (see §5.3 of IChad 

ioi^). 
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Figure 1. Velocity dispersion functions (VDFs) o f galaxies a t 
2 = and z = l. The VDF s at z = are t hose f r omlChael 1I2OIOI') 
and iBernardi et al.1 hOld ). The VDF by IChaj 1 I2OIOI) is a sum 
of the early-t ype and late-type VD Fs based on SDSS DR5 data. 
The VDF bv lBernardi et al.1 ||201(J) is a direct fit for all galaxies 
based on SDSS DR6. The VDF bv lCha3 toidti is also reproduced 
as a dashed curve on the upper right panel for comparison. The 
constraints on the VD Fs at 2 = 1 are based on strong lensing 
statistics described in IChael [2010'). Notice that strong lensing 
data probe only the range of 95km s^^ J; J; 300km (see 
the texts in § 2). 



Notice that the strong lensing surveys used to constrain 
the evolution of the VDF are limited to lensing galaxies with 
image splitting grea ter than 0.3 arcsec in the redshift range 
of 0.3 < z < 1 (see lChadlioioh . This lower limit on image 
splitting implies that the constrained evolution of the VDF 
is strictly valid only for a > 95 km s~^. This in turn cor- 
responds to > 1O^°-^M0 and Mvir > 10" '^Mq as wiU 
be shown below. Furthermore, although the surveys do not 
have physically meaningful upper limits on image separa- 
tions, the surveys (because of the small sample sizes) have 
only identified lensing galaxies with measured or implied 
stellar velocity dispersion a < 300 km s~^ (with correspond- 
ing < 1O" **M0 and Mvir < lO^^ '^Mo as shown below). 
This is why the constraints on the VDF evolution become 
weak toward large a as shown in Fig. [1] Hence any results 
from this work outside the above ranges must be regarded 
as extrapolations. 

For the HMF we use a typical numerical result from N- 
body simulations under the current ACDM cosmology while 
for the SMF we use the results from some representative 
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Figure 2. Loft panel: A typical mass func tion for dark haloe s as is 
produced from the ACDM simulation bv lReed et al ] l|2007l '). The 
adopted cosmological parameters are (OmOi^Ao) = (0.25,0.75) 
and as = 0.8. The displayed function has been corrected to in- 
clude subhaloes (see the texts in Section 3). Right panel: Observed 
ste llar mass functions . Solid curves are from the COSMOS survey 
by lllbert et al!* (2010') whi le dashed cu rves are from the Spitzer 
survey bv IPcrcz-Gonzalcz et ahl | |200^ . Notice that the Spitzer 
results show stellar-mass-downsizing evolution of the SMF while 
the COSMOS results do not. 



galaxy surveys. The HMF and the SMF are shown in Fig. [2] 
and more details are respectively given in Sections 3 and 4. 



3 CONNECTION BETWEEN THE 

OBSERVATIONAL VDF AND THE HMF 
FROM N-BODY SIMULATIONS 

In the ACDM hierarchical structure formation picture the 
dark halo mass function (HMF) evolves over cosmic time as 
a consequ ence of h ierarchical merging (e.g. IWhite fc Ree3 
[l978; L acev fc Colelll993 ). Accordingly, the statistical func- 
tions of galaxies such as the VDF and the SMF are also ex- 
pected to evolve. However, baryon physics complicates the 
evolutions of the VDF and the SMF making it non-trivial 
to compare the evolutions of the HMF, the SMF, and the 
VDF one another. Conversely, careful analyses of the coevo- 
lution of the HMF, the SMF, and the VDF may reveal key 
insights into galaxy formation and evolution processes. Here 
we compare the evolution of the VDF described in Section 2 
with the evolution of the HMF from cosmological N-body 
simulations. A comparison between the VDF and the SMF 
is given in the next section. 



Press & Schechteij Il974l: ISheth & Tormen 


1999,, 200I) or 


through A'^-body simulations (e.E. Ijenkins et al.1 


2001; 


SurinEel et al.1 l2005l: IWarren et all 20061: Reed et al.l 


2007; 


Lukic et al.l 20071: Tinker et al.l 20081). Recent high resolu- 



tion A^-body sim ulations have determined t he HMF reliably 
(e.g.lWarren et al 2006; Reed et al. 2007; Lukic et al.ll2007l ; 



[Tinker et alJbOOsh . For the HMF we g enerate 



a numeri- 



cal function using the code provided by iReed et al.l (|2007t ) 
taking the following cosmological parameters; f2mO = 0.25, 
Oao = 0.75 and erg = 0.8 consistent with the WMAP5 
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data l|Dunklev et al.ll2009l ). This function includes only dis- 
tinct haloes that are not parts of larger haloes. We correct 
it to include subhaloes since they may also host ealaxiesQ 
We use the simulation results by IConrov et alj (|2006 ') for 
the number fraction of subhaloes (/sub) as a function of 
maximum circular velocity. By relating the maximum cir- 
cul ar velocity to t he ha lo virial mass using the scaling given 
bv iKlvpin et Zl we find a varying fraction from 

/sub « 0.25 (r; 0. 18) at Myj^ < 10" M(7) (this mass scale 
corresponds to the lConrov et al" (|2006l ) completeness limit) 
to « 0.08 (« 0.08) at Mvir > IO^Mq for 2 = (z = 1). 
Notice that the mass of a distinct halo refers to the epoch 
under consideration while that of a subhalo is the mass at 
the time it accreted onto a larger halo. 

The mass of a halo (Mvir) may be linked to the stellar 
velocity dispersion (a) of the central galaxy for those haloes 
that host galaxies. If the halo did not host a galaxy in its 
centre, the central velocity dispersion (of dark matter par- 
ticles) would be entirely due to the dark mass potential. In 
reality, the central galaxy contributes to the central gravita- 
tional potential with the degree of contribution varying from 
one system to another. The functional relation a{Mvir) will 
depend not only on the stellar mass distribution of the re- 
siding galaxy but also how the dark halo has been modified 
due to the baryonic physics of galaxy formation. Moreover, 
the stellar mass distribution itself is correlated with Afvir 
to some degree. Thus, we may use the abundance matching 
relation between Mvir and a to gain new insights into the 
structure of the baryon-modified dark halo and the dynam- 
ical aspect of galaxy formation and evolution. 



3.1 The Mvir-cr relation at z = 

Fig. O shows the abundance matching M^ii-a relation at 
2: = 0. It shows both the relations ignoring intrinsic scatters 
and those taking into account an intrinsic scatter distribu- 
tion oi V = logjo((T/km s~^) as a function of Mvir. For the 
latter case the intrinsic scatter distribution is predicted by 
a bivariate distribution of a and A/* as a function of Mvir 
based on an observed scatter of logio(M*) at fixed Mvir and 
an observed scatter distribution of V as a function of M*. 
The reader is referred to Appendix A for a brief description 
and a following work (in preparation) for further details. 

The abundance matching relation is compared against 
the measured values of Mvir and a for individual galax- 
ies/clusters with z < 0.3. Although there are rmmerous 
galaxies/clusters for which either a or Mvir is reported, 
only for relatively few systems both Mvir and a have 
been measured reliably so far. First, we consider the best- 
studied Milky Way galaxy, for which r ecent m easurements 
appear to be reasonably concordan t ( Klvpin ct al. 20o3; 
Battaglia et al.ll2005l . l2006l : IXue et al.ll2008l ')PlThe data from 
Xue et all (2003) are displayed in Fig. [S] Second, we dis- 
play the results for 2 2 SLAGS lensing galax ies at mean red- 
shift of z ~ 0.2 bv iGavazzi et al.l (120071 ). They combine 
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Figure 3. The relation between the halo virial mass (Afvir) 
and the stellar velocity dispersion (cr) of the central galaxy at 
z = 0, inferred from the abundance matching of the velocity 
dispersion function and the ACDM halo mass function (see the 
texts in §3). The red curves are the results ignoring the scat- 
ter of V = logjQ((T/ km s"'^) at fixed Mvir- The black curves 
are the results taking into account the scatter shaded green. 
See Appendix A for a brief description of the scatter. The solid 
and d ashe d curves are respective ly based on the VDFs by IChad 
1I2OIOI) and lBernardi et al.l ll201(]|h The dotted line is a prediction 
of the SIS model of the halo. The curves are compared against 
the individual measurements for the following systems: Triangles 
- Milky Way (Xuc ct al. 2008 ). St ar - A weighted mean of 22 
SLAGS lenses l lGavazzi et al.l 120071) at 2 ~ 0.2 taking the stan- 
dard deviation of the individual mean values as the error on a. 
Circle - Lens sy stem Q0957+56 1 at z = 0.36. The velocity dis- 
persi on is from iTonrv fc Franxl |l99i) while the virial mass is 
from[Nakaj ima et al.| 1I2OO9I ). Squares - Galaxy cluster Abell 611 
l|Newma'i7et''an]2009l ) . The right square (gray) is based on the 
NFW halo model while the left square (solid) a generalised NFW 
halo model. 



strong and weak lensing and stellar kinematics to analyse 
the systems. Third, we consider the first ever discovered 
lens system Q0957-I-561 at z — 0.36, which is the best stud- 
ied lens system including a cluster for the lens potential. 
The velocity dispersio n for the central galaxy is reported by 
Tonrv fc Franx (119991) . The virial mass of the cluster is from 
Nakaiima et al.l (|2009l ) who derive the halo mass through 



^ This correction has only a relatively minor effect, in particular 
for large mass. 

^ The Andromeda galaxy is also a well-studied example, but the 
inferred virial masses appear to still vary by a factor of 2 (e.g.. 



iKIvpin et ai]|2002l : ISeigar et al]|200Sl) . 



weak len sing and find that the ir result is consistent with the 
result bv lGhartas et all (|2002l ) through X-ray observations. 
Finally, we consider g alaxy cluster Abell 61 1 that has been 
studied extensively bv lNewman et al.l (|2009l ) through a com- 
bination of strong and weak lensing and stellar kinematics. 
As shown in Fig. [3] these individual measurements are in 
excellent agreement with the Afvir-o" relation based on the 
abundance matching of statistical functions. This agreement 
bolsters the validity of the Afvir-"" relation at z = 0. 

Notice that the Mvir-o" relation shown in Fig. [3] is a 
curve rather than a straight line. Consequently, it does not 
match well the prediction by the singular isothermal sphere 
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Figure 4. The abundance matching A/vir-o" relation at 2 = 1 is 
compared with that at 2 = 0. It is consistent with zero evolution 
between z = I and z = for the strong lensing probed range 
o- < 300km s^i (A/vir < W^^-^Mq; see §2). 



(SIS) halo model (see, e.g., iLi fc Qstrike"3l2002l ). The failure 
of the SIS model is evident for Myir > 10^^ Mr as noticed by 
several authors (e.g. Li fc Ostrikeij|2002l : lKochanek fc Whitd 



l200ll : lBlumenthal et al.ll 19861 ). This result confirms the crit- 
ical halo mass Mc ~ 1O^^M0 below which the baryonic ef- 
fects start to become significant for the inner halo dynamics 
and structure. However, even for Mvir < IO'^Mq the SIS 
model is not very successful in matching the empirically de- 
termined (T(Mvir) curve. This implies that the SIS model is 
not precise as a 'global model' of the galactic halo despite 
the fact that a range of observational constraints support 
the isother mal profile for the inner part of the halo (see 
IChae|[2010l and references theirin). The underprediction of 
a by the SIS model for Afvir < lO^M© probably reflects 
the neglected concentration of the halo. The curvature in 
the Mvir-iJ relation may reflect the systematic variation of 
halo concentration but may also imply the varying bary- 
onic effects on the halo structures. The internal structures 
of the haloes may be constrained by combining dynamical 
constraints with the empirical Mvir-o" relation (in prepara- 
tion). 



3.2 The coevolution of the HMF and the VDF 

In Fig. 0] we compare the abundance matching Afvir-o" rela- 
tions at 2: = 1 and 0. This comparison shows little sign of 
evolution in the A/vir-o" relation for the strong lensing probed 
range a < 300km s~^ (Afvir < IO^ '^Mq; see §2). The near 
constancy in the Afvir-o" relation with z implies that the 
HMF and the VDF are coevolving in parallel. Namely, as 
the halo grows in mass over cosmic time, the central stellar 
velocity dispersion grows in accordance. The natural ques- 
tion to ask is then what the origin of this coevolution is. We 
discuss this in §5. 

Alternatively, we may transform the HMF into a VDF 
using the z = relation of Fig. |4] and assuming a certain 
evolution of the Mvir-o" relation. Fig.[S]shows the VDFs pre- 
dicted from the HMFs assuming zero evolution of the Mvir-o" 
relation. In Fig. [5] the HMF-converted VDFs are compared 
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Figure 5. The observationally constrained VDFs at z = 
and z = 1 are compared with the VDFs predicted from the 
ACDM halo mass function using the empirically determined rela- 
tion (T(A^vir) at z = assuming zero evolution (see Fig.|4}. There 
is a reasonably good match between the observationally derived 
VDF evolution and the halo-predicted VDF evolution. 



with the observationally derived local VDFs and the lens- 
ing constrained VDFs at 2: = 1. The VDFs at 2: = are in 
excellent agreement with each other. The VDFs at 2: = 1 
are also in agreement with each other. This exercise shows 
that under the simple assumption of the constancy of the 
Afvii-cr relation in time, the evolution of the HMF predicted 
by the current ACDM cosmology can match well the evolu- 
tion of the VDF constrained by strong lensing statistics for 
< 2 < 1. 



4 CONNECTION BETWEEN THE 

OBSERVATIONAL VDF AND THE SMF 
FROM GALAXY SURVEYS 

Many recent surveys of galaxies have been used to con- 
strain the evolution of galaxies through the LF or/and 
the SMF. The results are at variance. Many results ar- 
gue for relatively little evolution in the number den- 
sity of most massive galaxies and greater evolution of 
less massive galaxies over cosmic time, i.e. a "stel- 



Cimatti et all 12006 


; iFontana et all I2OO6I: iPozzetti et al. 


2007 


; Conselice et al. 20071; 


Scarlata et al. l2007l: Cool et al. 


200s 


; iPerez-Gonzalez et al. 


I2OO8I: iMarchesini et al.l |2009|). 



although there are results that do not particularl y sup- 

• Jb - 



port a mass-downsizing evolution (e.g., iBcll c t all I 2OO4I : 
iFaber et al]|2007l : lBrown et al.|[2008l : lllbert et aLlboiol ). The 

variance for the evolution of the SMF is not well understood 
but may be d ue to errors in measuremen ts and modelling of 
the SMF (see lLonghetti fc Saraccoll2009l ) and galaxy sample 
biases caused by cosmic variance (see, e.R. . lFaber et al ]|2007l : 
ICattaneo et al. I2OO8I : IStringer et aLlliooil for discussions). 

We have seen in the previous section that the strong 
lensing constrained VDF evolution is in line with the theo- 
retical HMF evolution. How well would the VDF match the 
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observed SMF from galaxy surveys? What would be the cor- 
relation between stellar mass (M*) and velocity dispersion 
(cr) and its evolution up to z = 1? The coevolution of the 
SMF and the VDF will depend on the evolution of the JVh- 
a relation. Hence the evolution of the VDF can be tested 
against the evolution of the SMF only when the M*-(t re- 
lation is known (or assumed) as a function of cosmic time. 
Conversely, by matching the observed SMF evolution from 
galaxy surveys with the VDF evolution from strong lensing 
statistics we may infer the evolution of the Ah-a relation. 

In the following we take the latte r appr oach. 

W e match the VDFs bv lChad (|201(]| 'l and lSernardi et all 
^2010) (their evolutions being constrained by strong lensing) 
with two total SMFs from galaxy surveys (see Fig. [2]) that 
are qualitatively different and are intended to encompass the 
current r ange of observations. One is the COSMOS SMF by 
lllbert et"a l. (2010) measured using 192,000 galaxies from the 
COSMOS 2-deg^ field. The COSMOS sample size is much 
larger than any other single data set th at has been u sed to 
derive the SMF up to z > 1. For the lllbert et al.] (|2010l ') 
SMF 2 = function is actually for 0.2 < z < 0.4 while 
z = 1 function is a mean of 0.8 < z < 1 . and 1.0 < z < 
1.2 functions. Notice that the lllbert et al.l HoT^) SMF does 
not show a stellar m ass-downsizing evolution. Th e other is 
the Spitzer SMF bv iPerez-Gonzalez et all (|2008l ')PI This is 
a ty pical SM F that shows a downsizing behaviour. For the 
IPerez-Gonzalez et all l|2008l ') SMF z = function is actually 
for < z < 0.2 while z = 1 function is a mean of 0.8 < 
z < 1.0 and 1.0 < z < 1 .3 functions. We note that the 
IPerez-Gonzalez et al.l (|2008| ) sample covers a sky area of only 
~ 664 arcmin^ and contains ~ 28, 000 sources for < 2 < 4. 



4.1 The M*-(T relation at 2 = 

Fig. [6] shows several examples of the abundance matching 
relation b etween and a at z = 0. For 2 = only we con- 
sider the iBernardi et~aLl l|2010l ') SMF as well as the COS- 
MOS and the Spitzer SMFs. These results have been ob- 
tained taking into account the effect of an adopted intrinsic 
scatter (the region shaded green) for V = log^Q{a / km s~^) 
of 0. 115 - 0.039 X (m - 1 0) with m = logi(,(Af*/M0) 
from iDesroches et al.l l|2007l ) (see Appendix A) . The details 
on the effect of the intrinsic scatter can be found in Ap- 
pendix A. Notice that the abundance matching M^-a re- 
lations have mild curvatures. For a linear approximation 
V = hm -\- const, the slope 6 varies from b — [0.23, 0.34] for 
m > 11.5 to 6 = [0.35,0.49] for m < 10.5. The abundance 
matching relations for all galaxies are compared with the di- 
rectly measur ed median relations for early-type galaxies in 
the literature (IDesroches et al.ll2007l : iHvde fc Bernardl2009l : 



IShankar etakl 120101 ). The abundance matching relations 
agree well with the early-type relations for M* > 10^'^ '^Mq. 
However, as decreases, the abundance matching rela- 
tions deviate systematically and increasingly from the early- 
type relations. This is expected and can be well understood 
by the fact that the late-type relation is different from the 



* iPerez-Gonzalez et al] ||2008| ') adopt the Salpeter IMF to calcu- 
late their stellar masses. To convert their stellar masses to those 
based on the Chabrier IMF, we divide by 1.7. 
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Figure 6. The abundance matching relation between the stellar 
mass (Mi,) and the stellar velocity dispersion (cr) of galaxies at 
2 = 0. Black (gray) solid, dashed, and dotted c u rves are respec- 
tively the results of matching the IChaei ||201C|) j Bernar di et al.l 
I2OI0I) VDF with t he CO SMOS SMF bv 'll bert et al.l 1 I2OIOI) , 
the Spit zer SMF bvlPerez'-G onzalcz ct al. (200^), and the SDSS 
SMF by IBernardi et al.l 1I2OIOI ) for all galaxies. These abundance 
matching results have been corrected for the effects of the in- 
trinsic scatter shaded green. The intrinsic scatter and its effects 
are described in Appendix A. The red solid, dashed and dot- 
ted curves/li nes are the measured m edian Mi,-a relat i ons re - 
spe ctively bv IDesroches et al.l ||2007| ). iHyde fc Bernardil {Sooi), 
and lShankar et al .1 l|20ld ) for early- type galaxies. Notice that the 
abundance matching relations agree well with the measured early- 
type relations at high A/* but deviate systematically as Af* gets 
lower because the late-type contribution becomes increasingly 
larger. 



early-type relation and the late-type contribution to the to- 
tal relation increases as Mi, decreases. 



4.2 The A'h-G relation at 2 = 1 and its evolution to 

2 = 

In Fig. [21 the abundance matching Mi,-a relation at 2 = 1 
is shown and compared with that at 2 = 0. The relation 
at 2 = 1 is also compared ag ainst the individual data 
poin ts for 0.7 < 2 < 1.3 fro m Idi Sereeo Alighieril (|2005l ) 
and Ivan der Wei et all (|2005l ). The resulting evolution of 
the Mi,-a relation varies depending mostly on the adopted 
SMF. Notice that the strictly valid range probed by the 
data is 1.97 < y[= logjQ((T/km s~^)] < 2.47 corresponding 
to 10.2 < m[= logic (A^^^/Mq)] ~ 11-8 (see §2). The results 
outside this range are extrapolations. 

For the COSMOS SMF the A/*-o- relation is consistent 
with zero evolution. The relation at 2 = 1 is also consistent 
with the measured data points. On the other hand, for the 
Spitzer SMF (a typical downsizing SMF) the 2 = 1 relation 
deviates systematically from the 2 = relation in particular 
for massive galaxies with Af* > IO^'^Mq. This appears to be 
the case for any downsizing SMF. This is because the evolu- 
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Figure 7. The M^^-a- relation at z = 1 is compared with 
that at 2; = 0, inferred from the abundance matching of the 
strong len sing constrained evol utions of the VDFs of IChaj 
l|201Clli and iBernardi et alj ||2010|) with the observed SMFs from 
galaxy surveys (see the texts in §4). T wo SMFs are used: the 
C OSMOS SMF by lllbert et al.l hout) and the Spitzer SMF 
by iPerez-Gonzalez et al. I | |2008D . The blue data points are the 
weighted means and their errors (thin error bars are dispersions) 
for mass intervals of 0.4 dex ba sed o n 47 galaxies for 0.7 < z < 1.3 
from ldi Serego Aliehieril bOOSi ') and lvan der Wei et al.l ||2005|') . All 
stellar masses are for the Chabrier IMF. 



tionary behaviour of a downsizing SMF is dissimilar to that 
of the VDF. This implies that a downsizing SMF requires 
a differential evolution in the M*-a" relation. In particular, 
according to the downsizing SMF a has to be lower at z = 1 
than z = at fixed (> 10"-^"" ''Mo). 

Then, which of the above cases (the non-evolving or the 
evolving case) of the K'h-o relation would be more consistent 
with other independent results on the structural evolutions 
of galaxies? 



4.2.1 Comparison with observed structural evolutions of 
galaxies 

According to recent studies on the structural evolu- 
tions of galaxies, there are observational indications 
that galaxy size evolves at fixed stell a r mass (e.g. 



Truiillo etal] l2007l: fvan der We i et all [200i: ICimatti et all 
20081 : Ivan Dokkum et al.. ,2008 : .Bezanson et all |2009| ) . How- 
ever, more recent studies find that physical mass den- 
sities (as opposed to effective densities) evolve little 
l|Hopkins et al.ll2009al : iBezanson et al.1 120091 ). This implies 
that velocity dispersion might evolve li ttle at fixed stellar 
mass. Indeed, ICenarro fc Truiilld l|2009l ) find a slow evolu- 
tion of (j(M*) since 2 ~ 1.6 from an analys i s of spectra 
of massive galaxies (see also ICappellari et al.|[2009l ). Nev- 



ertheless, ICenarro fc Truiilld l|2009l ) find a higher cr at a 
higher z for 0.5 x IO^^Mq < < 2 x 10"Mq (an evo- 
lution from a ~ 180 km s~^ at 2 ~ to ~ 240 km s~^ at 
1.6). However, as shown in Fig. [8] a similar data set actu- 
ally appears to indicate no evolution at all. The velocity 
dispersion at a fixed stellar mass of M* = 10^^ Mq rather 
than a range stays constant at cr ~ 210 km s~^ between 
2 ~ and ~ 1.8. Notice that we only use galaxies with 
10.75 < logiQ(M*/MQ)] < 11.25 and we convert the mea- 
sured value of cr at the measured value of for each galaxy 
to that at M* = IO^^Mq using an empirical relation found 
in Fig. [S] This prescription largely removes any systematic 
error arising from the differences in the stellar masses of 
the galaxies in the different redshift bin s. Another difference 
between the ICenarro fc Truiilld ll2009l) aiialysis and ours IS 
that for local galaxies ICenarro fc Truiilld (|2009l ) use SDSS 
DR6 data to derive stellar masses while we use only SLAGS 
galaxies for which two independent stellar mass mea sure- 
ments are available based on SDSS (|Grillo et al.ll2009l ) and 
HST (|Auger et al.ll2009^ photometric data. 

Incidentally, all the abundance matching results shown 
in Fig.[7]imply no or little evolution of a{z) at M* = lO^M© 
regardless of the VDF or the SMF used. Hence, our abun- 
dance matching results are in excellent agreement with, but 
at the same time are not distinguished by, our analysis of 
of the data from the literature for individual galaxies for 
^ ^ ^ 2. However, it is important to notice that none 
of the current observational results on the structural evo- 
lutions of galaxies indicate a negative evolution in cr with 
2 for any M*. Observational indications (e.g. stellar mass 
densities) are such that velocity dispersions cannot be lower 
at a higher 2 in case of evolution. Hence, according to our 
abundance matching results the downsizing SMF is i ncon- 
sistent with the VDF evolution from strong lensing bv lChad 
(2010) because it requires a lower o" at a higher 2 for massive 
galaxies. 

4.2.2 Comparison with predictions from cosmological 
hydrodynamic simulations 

At fixed stellar mass cosmological hydrodynamic simula- 
tions from the literature also find s low or little evolutions 
of velocity dispersion. For example, iHopkins et al.l l|2009bl ) 
combine dark halo merging with hydrodynamic simulation 
results and observed empirical properties of galaxies to find 
little e volutions of o-{z) at fixed Af* (Fig. IS]!. iHopkins et al.l 
find slow evolutions through more realistic cosmo- 
logical simulations taking into account various effects in- 
cluding equal and minor me rging, adiabatic expansion and 
observational effe cts (Fig.O. Cenarro fc Truiillol (|2009l ) pre- 
dict based on the lHopkins et al.l (l2009bl ) model a somewhat 
greater evolution using their analysis of spheroid size evo- 
lutions (Fig. |8}. These simulations are broadly consistent 
with the constraints from the current data as analysed above 
(Fig. [8} and the abundance matching results (Fig. [7]) in 
the sense that the predicted evolutions are slow and can 
be made in principle to ag ree with t he observational con- 
straints. Furthermore, Hopkins et al.l |2009b) find that the 
evolution of o-{z) at fixed has little sensitivity on M* for 
IO^Mq < M* ^ 1O^^M0. These simulation results are con- 
sistent with our abundance matching results based on the 
non-downsizing SMF but not with the downsizing SMF. 
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Figure 8. Upper panel: The observed stellar velocity dispersion 
(c) as a function of z at fixed stellar mass M* = IO^^Mq. The 
data points are based on galaxies with measured stellar masses 
in the range 10.75 < m[= logioCA^t/Mo)] < 11.25. Furthermore, 
to estimate the velocity dispersion at m = 11 as precisely as 
possible we use an empirical relation of V = bm + const where 
V = logj^g((T/km s~^) and we take b = 0.34 as estimated from 
Fig. [6] for 10.75 < m < 11.25 (the results are insensitive to the 
exact value of b for b > 0.2). The references for the d ata are as fol- 
lows: (1) z < 0.15, 8 (or 17 fo r bl ue point) galaxies — iBolton et al.l 
l|2008l 'l & iGrillo et al.l (Hool) (or lAuger et al.) (12003)') (2) 0.15 $ 



2 < 0.4, 9 (or 11 fo r blue point) gala xies -iBolton et all l|2008l ) 
iGrillo et al.l (120091') for lA ugcr ct ai] tOO^ )) (3) 0.6 < 2 < 0.8, 7 
g alaxies -Ivan der ^yel et al., (■2005i') (4) 0.8 ^ z < 1.0, 10 galaxies 
- Ivan der Wei et al.| ||2005^. Idi Sereeo Alighieril ll2005ll (5) 1.0 < 



z < 1 .3, 8 galaxies — Ivan der Wei et al.l (12001)1 di Se rego Alighieri! 
(I2OO5I ') (6) 1.6 ^ z < 1.8. 7 galaxi es -ICappellari et a l. (2009), (7) 
z=2.186, 1 galaxy (red point) - Ivan Dokkum et al.l (120091') . The 
solid line is the best-fit in the least-square fit of the data points 
and the dashed lines represent the errors in the slope. The single 
galaxy at z=2.186 is not used for the fit but consistent with the 
fit result at the 2cr level. Lower panel: The evolution factor of the 
velocity dispersion as a function of z. The black solid and dashed 
lines are the fit results fr om the upper panel. T he light hatched 
area is the prediction by iHopkins et al.l (l2009bl). Th e light gray 
area is th e result based on the iHopkinsetal .1 l|2009bl 'l model pre- 
sented bv lCenarro fc Truiillol ( |2009|) who use local SDSS data and 
similar data for 2 > 0.5 as used here but without converting the 
measured velocity dispersions to the values at a fixed stellar mass 
(they use a broad range 0.5 X 10" M© < Af^ < 2 X 10"M , 7)). T he 
dark gray solid curve is the prediction bv lHopkins et ahl (l2010h . 



To sum up, our abundance matching results based on 
the non-downsizing SMF are broadly consistent with obser- 
vational constraints on the structural evolutions of galax- 
ies and cosmological hydrodynamic simulation results. How- 
ever, the downsizing SMF gives a differential evolution in 
cj{K'h) with 2 that would not be consistent with observa- 
tional constraints or simulation results. 



5 IMPLICATIONS FOR THE ACDM 

PARADIGM AND GALAXY EVOLUTION 

At the heart of the current hierarchical structure forma- 
tion theory is the bottom-up build-up of dark matter haloes. 
Given that galaxies are believed to be born and centered in 
those haloes, what would be the evolutionary patterns of 
galaxies? Unlike dark haloes, galaxies have two distinctive 
properties, namely, the photometric property and the dy- 
namical property. Hence, there are two evolutionary proper- 
ties of galaxies to be considered. The star formation history 
of galaxies gives rise to above all the evolutionary patterns in 
the luminosity and stellar mass functions of galaxies. Most 
cosmological observations have been devoted to the photo- 
metric properties. Galaxy formation models, whether semi- 
analytical or hydrodynamical, have also tried to reproduce 
the photometric properties rich in observational data. No- 
tice that dark haloes do not have such photometric proper- 
ties. This means that connections between the photometric 
properties of galaxies and dark haloes are indirect and chal- 
lenging. Galaxy evolution in its dynamical property, which 
is the focus of this work, may be characterized by the evo- 
lutionary patterns in the velocity (dispersion) functions of 
galaxies. One may wonder whether the dynamical property 
of galaxies may be more intimately linked to dark haloes 
than the photometric properties do. How is the dynamical 
property of galaxy evolution related to dark haloes? What 
would be the role of the photometric property of galaxy evo- 
lution in this context? 

Current galaxy formation models cannot yet reliably 
predict central dynamical properties of galaxies. In this work 
we have compared the lensing-constrained evolution of the 
■VDF with the ACDM predicted evolution of the HMF and 
the observed evolution of the SMF. The main result is that 
the halo virial mass (Mvir), the galaxy stellar mass (M*) 
and the central line-of-sight stellar velocity dispersion (cr) 
are positively coevolving for the probed redshift range of 
< 2 < 1. 'What are the implications of the resuhs from this 
work for galaxy formation and evolution under the ACDM 
hierarchical paradigm? 

Halo mass-velocity dispersion relation: 'We find that 
Mvir-o" relation does not evolve between 2 = 1 and z = Q 
(Fig. [4]) for the entire probed range of halo mass from galac- 
tic haloes to cluster haloes. This empirical finding is insen- 
sitive to the choice of the HMF and the 'VDF from current 
simulations and observations under a concordance ACDM 
cosmological model. 

Pure dark halo simulations predict that a halo of mass 
Mvir at 2 = 1 is smaller (i.e. smaller -Rvir) but less concen- 
trated (i.e. smaller Cvir) than that at 2 = 0. As to the central 
velocity dispersion the two effects are opposite so that we 
may expect (j(Mvir) to evolve little as far as pure haloes 
are concerned. Let us consider this quantitatively using a 
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simple model. Without dissipational galaxy formation, the 
evolution of the central velocity dispersion would be primar- 
ily determined by the evolutions of the virial mass (Mvir), 
the virial radius (rvir) and the concentration (cvir)- The ve- 
locity dispersion is expected to increase (decrease, increase) 
if Mvir (fvir, Cvir) iucreases while the other two parameters 
are held constant. Cosmological A'^-body simulations predict 
that all three parameters (i.e. Afvir, ^vir and Cvir) increase as 
cosmic time evolves forward. Suppose a — v^ir f {cvir , r / r^ir) 
where a is the velocity dispersion in the central region (e.g. 
within O.Olrvir), v^ir = a/ GMvir/fvir is the circular velocity 
at the virial radius, and /(cvir, ?"/'"vir) is a model-dependent 
factor relating the two. A/^-body simulations show that Mvir, 
Tvir and Cvir all increase roughly b y a factor of 2 from 2 = 1 
to (see, e.g.. lWechsler et al.ll2002l ). Then, Vvir stays roughly 
constant and /(cvir, '"/f'vir) increases by about 15% from 
z = 1 to for an isotropic NFW model (see lLokas fc MamonI 
[2003). Hence we expect some enhancement in the velocity 
dispersion (i.e. a positive coevolution) in the course of the 
hierarchical growth of a pure dark halo from z = 1 to 0. 

For realistic haloes hosting (dissipationally formed) 
galaxies hydrodynamic simulations can be used to pre- 
dict the evolution of o"(Mvir). Unfortunately, current hy- 
drodynamic simulations do not pr edict robustly the bary- 



onic effects on halo structures (e . g. ["Blumenthal et all 



Gnedin et all '2004 : A badi et all l2009l : iTissera et al.l 
Feldmann et al., 2010.) 



198e 



201C 



Specifically, recent hydrodynamic 

simu l ations overpredict a at a given Mvir (e.g.ll issera et al.l 
I2OIOI : iFeldmann et aLllioioh . 

The finding that the Mvir-o" relation does not evolve 
for ^ 2: ^ 1 offers new insights into galaxy formation 
and evolution. It implies that the dynamical property of the 
central galaxy of a halo heis little to do with its history but 
is dictated by the final halo virial mass at least since z — 
1. Remarkably, this is the case for all haloes probed (with 
a > 100 km s~^). Implications of this finding are discussed 



consistent with zero evolution in a {Mi,) between z = 1 and 
0. On the other hand, the relation based on the Spitzer SMF 
(a typical downsizing SMF) implies a differential evolution 
in cr(M*): for M* > IO^^Mq the implied evolution in cr 
with redshift at fixed M* is negative while it is positive for 
M* < 10" M0. This means that based on the downsizing 
SMF a galaxy at 2 = 1 would have a shallower (steeper) 
mass profile than the local counterpart of the same stellar 
mass for KU > 10"Mq (A/* < 10"A'/q). 

How the above results on the evolution in (j(M*) are 
compared with other independent results on the structural 
evolutions of galaxies? First of all, we find little evolution in 
a{z) at M* = 10"Mo for < 2 < 1.8 from a careful analy- 
sis of the data in the literature (see Fig. [8|. This is in excel- 
lent agreement with the above abundance matching results. 
However, it cannot unfortunately distinguish the abundance 
matching results because M* = 10^^ Mq happens to be the 
critical mass for the downsizing SMF at which the evolution 
changes the sign. Second, many observational studies find a 
negative size evolution of galaxies with redshift implying a 
more steeply declini ng stellar mass profile at a higher 2 (e.g. 
| Trui il lo ct al. 2007,: [van der Wei et al.1 12OO8I : ICimatti et all 
lioosi :! van Dokkum et al.ll2008l ') . However, more recent stud- 
ies find that stellar mass density profiles of the inner regions 
up to several kilo-parsecs are consisten t with no evolution for 
massive galaxies with A/, > 10^^ M© (jHopkins et al.|[2009al : 
[Bezanson et al. 2009 ) . According to these studies, however, 
it is not clear whether stellar mass density profiles evolve 
beyond the inner regions. Whatever the case these results 
can only imply a similar or larger a at fixed M* at a higher 
redshift contradicting the abundance matching results based 
on the downsizing SMF. 

What do hydrodynamic simulations pre dict on the evo- 
lution of the relation between Af* and a? iHopk ins et al.l 
(2009b) combine galaxy merging with hydrodynamic sim- 
ulation to find a little evolution of cr with 2 at fixed M* for 



below in the context of the coevolution of Mvir, cr and Af*. any 10^ M a < M* ^ lO^^M©. In particular, iHopkins et al.1 



Stellar mass-velocity dispersion relation: The A'h-a re- 
lation at 2 = shows a power-law relation M* oc a''^" 
with a varying power-law index 7sm ranging from [2.9, 4.4] 
for A'/* > 10^^-^Mq to [2.0,2.9] for M* < 10^°-^Mq. Let us 
compare the Mi,-a relation with power-law correlations be- 
tween luminosity and internal velocity parameter, namely 
the TuUy-Fisher relation for the late-type population and 
the Faber- Jackson relation for the early-type population. 
The observed TuUy-Fisher r elation exponent 7t f lies be- 
tween 2.5 and 3.5 (see §2.3 or lPizagno et al.ll2007l '). The tra- 
ditional value for the Faber- Jackson exponent 7fj for early- 
type galaxies is ~ 4. However, an extensive analysis of SDSS 
DR5 early-type galaxies reveals that 7fj varies systemati- 
cally fro m 2.7 ±0.2 at Lt to 4.6 ± 0. 4 at the upper luminos - 
ity end (|Choi et al.1 l2007l : see also iDesroches et al.1 120071 ). 
The abundance matching M*-cr relation for all galaxies 
can match well these Fab er- J ackson / TuUy-Fishe r relations 
in co njunction with measured Af*/L ratios (e.g. iBell et al] 
120031 ). 

In Fig. [7] the A'U-a relation at 2 = 1 is compared with 
that at 2 = based on two VDFs and two SMFs that are 
meant to encompass the current range of observations. As 
can be seen in the figure, the implied evolution depends 
sensitively on the adopted SMF and to a less degree on the 
adopted VDF. The relation based on the COSMOS SMF is 



2009bl ) explicitly predict t hat the M^-g re l ation evolves lit- 
tle between 2 = 1 and 0. i Hopkins et al] (|2010l ) take into 
account a number of possible effects in their cosmological 
simulations and find slow evolutions of a with 2. These re- 
sults are consistent with the evolution in cr(M*) with 2 based 
on the COSMOS SMF but not with the downsizing SMF. 

VDF evolution: concord or conflict with observed galaxy 
evolutions? We have already compared the evolving VDF 
with the evolving HMF and the evolving SMF. We find that 
the evolutions of the HMF, the VDF and the SMF appear 
concordant, but that the downsizing SMF is disfavoured 
because its implied structural evolutions are unlikely. It is 
clearly worthwhile to put the VDF evolution in a broader 
context of recent cosmological observations on galaxy evo- 
lutions. 

The lensing constrained VDF evolutions show that 
the number density of massive early-type galaxies (cr > 
220 km s~^) not only evolves significantly but also shows 
a differential evolution (see Fig. [l}: the higher the ve- 
locity dispersion, the faster the number density evolu- 
tion (the "velocity-upsizing" behaviour), probably mean- 
ing an "mass-upsizing" behaviour (i.e. the behaviour of 
more massive galaxi e s asse mbling later in cosmic time). 
iMatsuoka fc Kawaral (|2010l ) has just recently compiled a 
large number (~ 60, 000) of massive galaxies (M* > 
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lO^M©) over a large sky area (55.2 deg^) from the UKIRT 
Infrared Dee p Sky Survey (UKIDSS) and the SDSS II Super- 
nova Survev. lMatsuoka fc Kawaral l|201(t l find a significantly 
greater number density evolution for M* > than 
M* < fO^^'^M0 out to z = 1 consistent with the hierarchical 
evolution. The parallel evolution of the VDF and the SMF 
would imply no evolution in the tot al mass profile of galax- 
ies. Indeed, strong lens mo delling ( Koopmans et all 1200' 
ISaha et al ] |2007l : IWinn et al.. 2004 : Treu fc Koopmanslf20& 
and velocity dispersion measurements to a high redshift 
(Fig. |8} support non-evolution in total mass profiles even 
if stellar mass profiles evolve. Galaxy merging is an inde- 
pendent route to probe the build-up of galaxies over cosmic 
time ( e.g.lWhite et aL|l2007l : iMasiedi et aLl200^ : IWake et all 



l2008l : Ide Ravel et all [2OO9I : Bundv et all f2009f i. The mere 
fact that thousands of merging events have been observed 
is the evidence for some sort of hierarchical mass assembly 
going on. The issue is the merging rate and its dependence 
on galaxy mass. Observed merging rates are at variance and 
cannot test the hierarchical assembly of massive galaxies ro- 
bu stly. It is, h owever, worth noting the more recent result 
bv lBundv et al . (2009) that merging rate is greater for mas- 
sive galaxies with M<- > 10^^ Mq than less massive galaxies. 
The lBundv et al.l (|2009D result is in line with the hierarchi- 
cal mass assembly and agrees qualitatively with the VDF 
evolution. 

Coevolution of A/vir, o" and M* and implications for 
the ACDM paradigm: In the above discussions we have con- 
sidered the connections of a with Mvir and M* separately. 
The results that a is coevolving in parallel with both Mvir 
and M* necessarily imply a similar coevolution of Mvir and 
Mi,. Fig. [9] shows the abundance matching Mvir-M* rela- 
tions at 2 = and 2 = 1. The results based on the COS- 
MOS SMF give little evolution in the Mvir-M* relation for 
Mvir ^ IO^^Mq implying a parallel coevolution of Mvir and 
Mi, with cosmic time. This is consistent with the little evo- 
lutions in the Mvir-cr and the M^-a relations based on the 
same SMF. Hence, we are left with the simple picture that 
Mvir, Af* and a are all coevolving so that a halo of given 
mass > 10^^ M0 has on average the same stellar mass and 
the same stellar velocity dispersion for its central galaxy in- 
dependent of redshift for < z < 1. Some indications of 
evolution in the Mvir-M* relation (and possibly in the Mi,-a 
relation) for Mvir < 10^^ Mq may imply a differential evolu- 
tion among Mvir, M* and a. We cannot address this issue 
for low-mass haloes here because the strong lensing con- 
strained VDF evolution from this work breaks down at low 
a. The results based on the (downsizing) Spitzer SMF give 
a mild differential evolution in the Mvir-M* relation at large 
massed as it is the case for M*-cr relation. However, in the 
above we have argued that a downsizing SMF is unlikely. 

How the parent dark halo is dynamically related to the 
residing galaxy as a function of cosmic time is a fundamen- 
tal question for galaxy formation and evolution. In this work 
we have studied the connection of the halo mass (Afvir) with 
the stellar velocity dispersion (a) and the stellar mass (M*) 
of the central galaxy for < 2 < 1. According to our results 



^ Given the observational uncertainty of the Spitzer SMF at 
z = 1 it is only mar ginally inconsistent with zero evolution (see 
iBehroozi et al.ll2010l) . 



the stellar dynamics in the galaxy (characterized by a) is 
closely linked to the parent halo mass Mvir independent of 
redshift for < 2 < 1. A link between the central parti- 
cle velocity dispersion and Mvir is expected in the ACDM 
paradigm because it predicts on aver age a universal den- 
sity p rofile independent of Mvir (e.g. iNavarro et al.l |2004| . 
I2OI0I and references theirin) along with well-defined scal- 
ing relations of th e structural parameters with Mvir (e.i 



iBuUock et al]|200ll : iMaccio et aLll2007l : iKlvpin et al]|201ol ). 

According to the coevolution, dissipational baryonic physics 
involving star formation that may have modified the central 
potential of the halo has not destroyed but appears to have 
refined the link. 

It is then natural to suggest that dissipational baryonic 
physics results in on average a rescaled universal (or univer- 
sal class) density profile of the stellar plus dark mass distri- 
bution, or perhaps more realistically a universal (class) den- 
sity profile of dark matter combined with a well-correlated 
class of stellar mass distribution. What would then the 
baryon-modified universal (class) density profile look like? 
The observation that the inner density profile of the g alaxy 
plus halo system is on average close to isothermal (see iChad 
l2010l for a review and references) combined with the expecta- 
tion that baryonic effects a re not likely to be important well 
outside the scale radii (e .g. iGnedin et al.|[20o3 : lAbadi et al.l 
'2009': 'Tis sera et al.ll2010l ) leads us to suggest a modified or 
generalized NFW (GNFW) profile in which the inner total 
density profile is close to isothermal (with a possible system- 
atic variation with Mvir) while the NFW profile is kept at 
large radii well outside the stellar mass distributions. The 
GNFW profile is then (in an average sense) preserved in the 
successive merging of GNFWs. Furthermore, the evolution 
of the concentration of such an GNFW profile with 2 at 
fixed Mvir conspires with the evolution of the virial radius 
rvir with z to lead to a non-varying a with 2. It is not well 
understood at present whether this is just a coincidence or a 
revelation of a fundamental mechanism in galaxy formation 
and evolution. It is also not clear whether the non-evolving 
Mvir — cr relation extends to a higher redshift, i.e., since when 
the HMF and the VDF have been coevolving in parallel in 
cosmic history. 

The observation that the amount of star forma- 
tion (i.e. stellar mass M i, ) is correlated with Mvir (e.g. 



'Conroy fc Wechsleil 
i2009; Behroozi et all 



200E 



Moster et all l2009l : ICuo et all 
2010|) is consistent with the above pic- 



ture. Namely, a larger halo undergoes a larger amount of 
star formation needed to modify the greater potential well. 
If Mi, were perfectly correlated with a at fixed Mvir, the 
correlation between Mvir and a would be just a by-product 
of the Mvir-M* correlation. However, although there is some 
good correlation between M* and a for all galaxies (i.e. re- 
gardless of their haloes), the correlation between M* and cr 
at fixed Afvir is weaker (in preparation). Hence, under the 
above picture the Mvir — cr correlation is originated from 
the ACDM haloes and the amount of star formation set by 
Afvir rescales the correlation. We then expect some correla- 
tion between M* and a at fixed Mvir because the boost of 
a depends on the degree of the baryonic effects on the halo 
characterized by M* (in preparation). 

Let us compare the coevolution and the above picture 
motivated by it with some pictures (or interpretations) and 
numerical simulation results of the ACDM paradigm that 
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LOGio(M,ir /Mo) 

Figure 9. The abundance matching relations between M^ii (halo 
virial mass) and M* (central galaxy stellar mass) at z = (black 
curves) and z = 1 (red curves) based on the COSMOS (solid 
curves) and the Spitzer (dashed curves) SMFs. These results are 
based on a constant intrinsic scatter of 0.16 for logj^g(M*/M0) 
at fixed A^vir- 



have been discussed in the literature. First, the continual 
growth of the central stellar velocity dispersion and the 
stellar mass accompanying the growth of the halo over cos- 
mic time would be inconsistent with a strictly "stable core 
concept" f or m assive galaxies (e.g. iLoeb fc Peebles! l2003l : 
iGao et al.l liooi ) even after 2 = 1. However, the growth 
slopes for a and M* are shallower for more massive haloes 
according to the abundance matching results (see Fig. 3] 
and Fig. [9]) . Hence a weakly evolving core of massive haloes 
would be consistent with our results. Second, the picture 
shares the concept of "u niversal density profile" with the 
attractor hypothesis (e.g. iLoeb fc Peeblej l2003l : iGao et al] 
[2004 ) . However, there is a clear distinction between the two. 
The attractor hypothesis assumes that the universal NFW 
profile is preserved or restored in hierarchical merging of 
haloes (hosting galaxies) while the present picture assumes 
that the baryon-modified universal total density profile (i.e. 
the GNFW profile) is preserved once it is created. The lat- 
ter property is supported by dissipation l ess merging simu- 
lations (e.g. iB ovlan- Kolchin fc Mai l2004l : iKazantzidis et al] 
l2006l : lNipotiet al. 200^^ 

Test of the ACDM paradigm?: The basic tenet of the 
ACDM structure formation theory is the hierarchical mass 
assembly. What the theory predicts is the distribution of 
dark matter haloes. Connecting observed galaxies with the- 
oretical dark haloes is a major goal of cosmological research. 
The difficulty of testing the ACDM paradigm arises from 
the complex physics of galaxy formation within the halo 
and the induced modification of the halo structure. A nec- 
essary condition for a successful model is to reproduce the 
basic statistical properties of the observed local Universe, 
such as the luminosity, stellar mass and velocity functions 
of ga laxies and their correlations (see iTruiillo-Gomez et all 
l201(lh . However, a successful reproduction of the z = sta- 
tistical properties of galaxies is not sufficient. A successful 



model must predict correctly the evolution of the galaxy 
properties. Semi-analytic models of galaxy formation have 
paid much attention on the galaxy luminosity and stellar 
mass functions. The current generation of these models can 
reproduce the z = functions reasonably well, but fail to 
match their observed evolutions (see, e.g.. iFontanot et al.l 
I2OO9I : IStringer et al]|2009l : ICattaneo et al lboosi ). 

The galaxy luminosity and stellar mass functions have 
much to do with the complex baryonic physics of star for- 
mations, AGN activities, feedbacks, etc. Hence the evolu- 
tions of the galaxy luminosity and stellar mass functions 
can only provide indirect tests of the underlying ACDM 
paradigm. The stellar velocity or velocity dispersion of the 
galaxy residing in the centre of a halo probes the grav- 
itational potential of the baryon plus dark matter sys- 
tem. Since the dark halo is expected to be modified in 
the course of the dissipational galax y forma t ion process 
(e.g.lBlumenthal et al."l986': ' Cnedin et al.i2004l : iRudd et al.1 
l200alAbadi et al.ll2009. : ,Tissera et al.ll2010l ) and the central 
potential is likely to be dominated by the baryonic matter, 
the velocity (dispersion) function evolution itself is not a 
direct probe of the ACDM paradigm either. However, the 
velocity (dispersion) function is separated from much of the 
baryonic physics but has only to do with its dynamical ef- 
fect. Hence, once the dynamical effect of galaxy formation 
is well accounted for, the evolution of the velocity (dis- 
persion) function offers an useful complementary probe of 
the structure formation theory. While it is challenging to 
measure reliably the evolution of the velocity (dispersion) 
function through conventional galaxy surveys, strong lens- 
ing statistics in a well-defined survey provides an excellent 
opportunity to constrain the evolution of the velocity (dis- 
pers ion) functio n through the image splitting distributions 
(see lChaellioiol ). Current strong lensing statistics is limited 
by the small sample size. However, future cosmological sur- 
veys including (but not limited to) the Dark Energy Sur- 
vey, the Large Synoptic Survey Telescope and the Square 
Kilometre Array will increase dramaticall y the number of 
strong lenses (see lOguri fc Marshall I2OI0I ) allowing to put 
tight constraints on the evolution of velocity (dispersion) 
functions. 



6 CONCLUSIONS 

Through an abundance matching analysis of the lensing con- 
strained VDF evolution along with the theoretical HMF and 
the observed SMF from galaxy surveys, we find the follow- 
ing. 

(i) The dark halo virial mass-central stellar velocity dis- 
persion (A/vir-o") relation at 2 = is in excellent agree- 
ment with the observed properties of low-redshift individual 
haloes. 

(ii) The stellar mass-central stellar velocity dispersion 
(M^-cr) relation at 2 = is consistent with the local scaling 
relations of galaxies in the literature. 

(iii) The Mvir-o" relation does not evolve between 2 = 1 
and independent of current observation and simulation 
data. 

(iv) The Mi,-a relation does not evolve between 2 = 1 
and for the COSMOS SMF. This is well in line with the 
observed non-evolution of a with 2 at = lO^^M©. This 
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is also consistent with the predicted little or mild evolution 
of a with z insensitive to M* from cosmological simulations. 
However, the Spitzer SMF (a typical downsizing SMF) re- 
quires the Mi,-a relation to evolve in a differential way that 
is not supported by independent observational results on the 
structural evolutions of galaxies in the literature. 

(v) The non-evolution in the Myir-o" and the M*-(j rela- 
tions imply a parallel coevolution of Mvir, and a between 
z = 1 and 0. This is corroborated by the little evolution in 
the abundance matching AI^ii-Mi, relation between z — 1 
and for Mvir > IO^^Mq. 

(vi) The parallel coevolution of Afvir, cr and Af* with z 
may imply a universality and regularity in galaxy formation 
and evolution despite complex baryonic physics processes. 
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APPENDIX A: INTRINSIC SCATTERS AND 
BIAS CORRECTIONS FOR THE ABUNDANCE 
MATCHING RELATIONS 

This work is concerned with the abundance matching (AM) 
relations between a and Mvir and between a and A'h. The 
AM is intended to recover the true median relation. If there 
were not any intrinsic scatter in the relation between two 
parameters, the AM by equation © of two statistical func- 
tions would recover the true relation exactly. In reality the 
distribution of two observables in a plane has intrinsic scat- 
ters around the median relation. When there are such in- 
trinsic scatters, the AM by equation ^ will give a biased 
median relation that is different from the true median rela- 
tion. Here we estimate the bias and correct the AM relation 
by equation ^ to obtain the corrected relation. The bias- 
corrected AM relation is then checked for self-consistency 
through a Monte-Carlo simulation. In other words, we esti- 
mate the bias so that the corrected relation reproduces from 
one statistical function to the other through a Monte-Garlo 
simulation based on the intrinsic scatter. 

For our purpose a knowledge of the intrinsic scatter 
is required. There has been no measurement or simula- 
tion for the intrinsic scatter in the Mvir-cr relation. On 
the other hand, there have been measurements for the 

Desroches et al.l l2007l : 



relation between M<- and cr (e , , . . 

iHvde fc Bernardil l2009l : IShankar et all |201G| ). Hence we 
study first the M*-cr relation using a model intrinsic scat- 
ter motivated from observed intrinsic scatters. We then give 
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m=LOGjo(MyMQ) F=LOGj(,(CT/km s ') 

Fi gure Al. Lef t pane l: The SMF for all galaxies measured 
by iBernardi et al based on SDSS DR6 and data points 

realised by a Monte-Carlo simu lation. Right panel: T he black 
curve is the VDF measured by IBernardi et al I l|201(]| ) for a > 
125 km s~^ based on SDSS DR6. The data points and green curve 
are the results of a Monte-Carlo simulation based on the bias- 
corrected AM relation and the intrinsic scatter shown in Fig. IA2I 
The red curve is the simulation result based on the biased AM 
relation shown in Fig. IA2I 



some results on the Afvii-cr relation that come from a pro- 
cedure to simultaneously derive an intrinsic scatter and the 
bias-corrected Afvir-o" relation by considering a bivariate dis- 
tribution of Mi, and ct as a function of Mvir (in preparation). 
Studies on the effects of the intrinsic scatters in the A/vir- 
M^, and the M^-^-L relati ons can be found respe ctively in 
iBehroozi et~al] (|2010l ') and IXasitsiomi et all l|2004 ). 

For the p urpose of demon s tratio n we use the SDSS SMF 
and VDF by IBernardi etHI l|2010l ) that are displayed in 
Fig. lAll For the intrinsic scatter of y[= logj^Q(cr/km s"'^)] as 
a function of m[= logj^Q(Af*/M0)], we adopt a linear model 
given by 0.115 — 0.039 x (m — 10) that is derived from the 
iDesroches et al.l l|2007l ) measurements of early-type galax- 
ies for_m^_>_J_0^^_and are co nsist ent with the measurem ents 
by iHvde fc Bernardil (|2009l ) and IShankar et al.l l|2010l ). For 
all galaxies including late-type galaxies, the intrinsic scat- 
ter will be more complicated than this. In this work we do 
not attempt to consider an intrinsic scatter distribution for 
all galaxies for the following two reasons. First, there have 
not been any published measurement results of the intrin- 
sic scatter for all galaxies. Second, the intrinsic scatter for 
early-type galaxies will match that for all galaxies at large 
stellar masses where the bias in AM is most significant. In 
other words, we can reliably estimate the greatest bias in 
AM using only the intrinsic scatter for early- type galaxies. 

Fig. IA2I shows the biased (red curve) and the bias- 
corrected (green curve) AM relations. The dashed curves 
around the bias-corrected AM relation represent the adopted 
intrinsic scatter described above. Fig. lAll shows the input 
VDF and the Monte-Carlo simulated VDFs from the in- 
put SMF based on the biased and bias-corrected AM re- 
lations and the adopted intrinsic scatter. Notice that for 
the bias-corrected AM relation the simulated VDF closely 
matches the input VDF. Hence the required self-consistency 
is gained. If we use different input SMFs and VDFs, we will 




Figure A2. Upper panel: Red curve is the AM relation given 
by equation [2] for the SMF and the VDF shown in Fig. lAll while 
green curve is an adjusted relation. Green dashed curves rep- 
resent la dispersion assuming a gaussian distribution of V[= 
logj^Q(o'/ km s~^)]. Data points have been realized from the SMF 
using the adjusted relation and the adopted scatter. A VDF de- 
rived from these simulated data points matches nearly perfectly 
the input VDF. The adjusted relation is referred to as the bias- 
corrected relation in the texts. Lower panel: The difference be- 
tween the initial AM relation (equation [2]l and the adjusted rela- 
tion, referred to as the AM bias. Red curve is that for the upper 
panel. Other curves are for the following input SMFs and VDFs: 
(1) black solid - Chae VDF/COSMOS SMF; (2) black dashed - 
Chae VDF/Spitzer SMF; (3) black dotted - Chae VDF/Bernardi 
SMF; (4) gray solid - Bernard! VDF/COSMOS SMF; (5) gray 
dashed - Bernardi VDF/Spitzer SMF; (6) gray dotted - Bernardi 
VDF/Bernardi SMF (identical to the red curve). 



of course get different AM biases. Several biases required for 
the z = Q SMFs and VDFs used in this work can be found 
in the bottom panel of Fig. IA2I To estimate the biases re- 
quired for the z = 1 functions we use the predicted VDFs at 
z = 1 based on the best-fit evolution parameters from strong 
lensing statistics along with the observed z = 1 SMFs. The 
magnitude of the biases for z = 1 is about the same as that 
for 2 = and we do not display the 2 = 1 biases. 

The same procedure would be followed for the Afvir-cr 
relation if there were an observational intrinsic scatter as 
for the Mi,-a relation above. Without a knowledge of the in- 
trinsic scatter for the Mvir-o" relation we devise a procedure 
that allows us to derive simultaneously the intrinsic scatter 
and the median relation for Mvir and a (in preparation). 
A full description of the procedure is beyond the scope of 
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this paper and the reader is referred to a following paper 
in preparation. Here we only give a brief description of the 
procedure and quote a simple result. The idea is to use a 
bivariate distribution of M* and a at fixed Mvir noticing 
that M* and a are expected to be correlated. Then, we de- 
termine simultaneously the scatter of a and the correlation 
coefficient between Mi, and a given the observed scatter of 
at fixed Afvir so that the resulting Mi,-G relation is con- 
sistent with the given relation based on observations. The 
results depend on the input scatter of M* and the input M*- 
o" relation. Fig. [4] shows a simple result based on a standard 
deviation of 0.16 for logj^Q(M*) at fix ed Myir and the above 
described M*-(t relation based on the lBernardi et all (|2010l ') 
observational results. 



APPENDIX B: 

Z = 



THE «v 



RELATION AT 



button et aD (|201(]| ') constrain the relation between Wvir and 
fopt through combining observationally derived M*-Mvir 
and Af*-«opt relations (nvir and «opt are the circular rotation 
velocities at the virial and the optical radii respectively). 
The abundance matching A/vir-o" relation from this work 
may be transformed into a Uvir-^'opt relation using an em- 
pirical relation between iiopt and a. The virial velocity «vir 
is defined b y -^/GA/vir/ryir where th e virial radius at z = 
is given by l|Brvan fc Normanlll998h 



Mvir 



10i2ft-iM, 



1/3 



kpc 



(Bl) 



for the adopted cosmology and we take h — 0.7. 

For early-type (elliptical and lenticular) galaxies, the 
direct estimate of Vo pt /cr ranges from » ^2 to ^ 1.7 (e.g. 
Courteau et~aIll2007l :lHo 2007.; Pizzella et al.ll2005l : lFerraresel 



20021 ). For late^tvpe galaxies the lChael (|2oTol ) late-type VDF 



has actually been transformed from a circular velocity func- 
tio n (VF ) assu ming Wopt/o" = \/2. Hence we can transform 
the lChael |2010''l late-ty pe VDF back to the original VF using 
the same factor. The Bernardi et al.l (fioiQ ') VDF is a directly 
measured function based on SDSS spectroscopy. It is impor- 
tant to note that the SDSS measured velocity dispersions 
for late-type galaxies with small bulges (or without bulges) 
come mostly from rotational motions (M. Bernardi, private 
communications), meaning that the se small-bulge ( bulge - 
less) systems are not missed in the iBernardi et al] (|2010l ') 
VDF. Hence we need an independe nt knowledge of Vo^t/a 
for SDSS galaxies to transform the IBernardi et al.l (|2010l ) 
VDF to a VF. Without a measured value of Wopt/o" for SDSS 
galaxies we must resort to other measurements. From the 
literature we find Vopt/<y ~ 1.4 — 2 for la te-type galaxies 
depending on the bulgc-to-disk ratio (e.g. 'Courte au et al.l 
[2007; Ho 2007; Pizzella et al. 2005; Fcrrarcsc 2002). 

Based on these literature values of fopt/o" for early- and 
late- type galaxies we adopt a range of «opt/o" = \/2 — 1.7 
independent of galaxy type to estimate the circular velocity 
function of galaxies and then the t^vir-t'opt relation through 
abundance matching. Fig. IBll shows the likely range of 
the median value for fopt/fvir as a functi on of Vvh. Com- 
pared with Fig. 5 of iDutton et al.l (|2010l ) our results for 
ma ssive galaxies ( tiyjr x *. 10^'^ « 200 km s~^) overlap with 
the lDutton et all (|2010l ') results for early- type galaxies. For 
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Figure Bl. The relation between the halo circular velocity at 
the virial radius (I'vir) and the stellar circular velocity in the op- 
tical region (iiopt) of the central galaxy at 2 = 0, inferred from 
the abundance matching relation between Myir (halo virial mass) 
and (T (central stellar velocity dispersion) from the VDF and the 
ACDM halo mass function (see the texts in §3.2). The solid and 
das hed curves are respecti vely based on the VDFs bv lChaj l|2010l ) 
and IBernardi et al.l 1I2OIC1I ) . For each set of the results the lower 
and upper curves are respectively based on t^opt = \/2a and 

''opt = 1.7(7. 



Vviv ^ 10^ '^ km s~^ our results give a scaling of t^opt oc u^j^. 
with 7 « 0.3 — 0.4. This scaling is consistent with previous 
results for massive galaxies based on strong lensing statis- 
tics and other methods (see lChae et al.|[2006l and references 
theirin) . 

However, our results for less massive galaxies {vvii < 
10^'^kms^^) give systematic ally higher median v alues of 
«opt/t'yir compared with the iDutton et al.l (|2010l ) results 
for late-type galaxies although the estimated intrinsic scat- 
ters (not shown here) overlap. Our results are consistent 
with declining rotation curves near th e viria l radii for late- 
type galaxies while the IDutton et al.l (|2010t ) results imply 
fiat rotation curves right up to the virial radii. The pos- 
sible discrepancy may imply some unidentified systematic 
errors in one or both of the results. For our results the 
possible sources of systematic errors include the adopted 
VDFs and the adopted relatio n between a and fop t- In or- 
der to be consistent with the iDutton et all (|2010l ) results 
the number densities of galaxies would have to be signif- 
icantly lower at relatively low velocity disp ersions. We use 
two independently d etermined VDFs (i.e. the lCha j|2010l and 
IBernardi et al.|[2010l VDFs) a nd they give similar re sults for 
Wopt/fyir. In fact, the adopted IBernardi et al.l (|2010l ) VDF is 
the modified Schechter fit result for a > 125 km s"'^. This 
VDF gives an und erestimate of ga l axy n umber densities at 
low a according to IBernardi et all (|2010l ) data. Hence if we 
were using th e raw galaxy numbe r densities, the discrep- 
ancy with the lDutton et all l|2010l ) results would get worse 
at low a. The adopted range Wopt/o" = V2~ 1.7 is suggested 
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by a broad range of observations. A possible cause of error 
may be that bulgeless late- type galaxies have Uopt/o" > 1.7 
so that we need to adjust the overall value higher. How- 
ever, if we adopted Vopt/o > 1.7, the discrepancy with the 
iDutton et al.l l|2010l) results would get worse. 



